Magnetic field— induced Kondo effects in Coulomb blockade systems 



M. Pustilnik", L. I. Glazman", D. H. Cobden'', and L. P. Kouwenhoven c 

a Theoretical Physics Institute, University of Minnesota, 116 Church St. SE, Minneapolis, MN 55455 
b Department of Physics, University of Warwick, Coventry, CV4 7AL, UK 
c Department of Applied Physics and ERATO Mesoscopic Correlation Project, 
Delft University of Technology, P. O. Box 5046, 2600 GA Delft, the Netherlands 

(February 1, 2008) 

We review the peculiarities of transport through a quantum dot caused by the spin transition 
in its ground state. Such transitions can be induced by a magnetic field. Tunneling of electrons 
between the dot and leads mixes the states belonging to the ground state manifold of the dot. Unlike 
the conventional Kondo effect, this mixing, which occurs only at the singlet-triplet transition point, 
involves both the orbital and spin degrees of freedom of the electrons. We present theoretical and 
experimental results that demonstrate the enhancement of the conductance through the dot at the 
transition point. 
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I. INTRODUCTION 

Quantum dot devices provide a well-controlled object 
for studying quantum many-body physics. In many re- 
spects, such a device resembles an atom imbedded into 
a Fermi sea of itinerant electrons. These electrons are 
provided by the leads attached to the dot. The orbital 
mixing in the case of quantum dot corresponds to the 
electron tunneling through the junctions connecting the 
dot with leads. Voltage V g applied to a gate - an elec- 
trode coupled to the dot capacitively - allows one to con- 
trol the number of electons N on the dot. Almost at any 
gate voltage an electron must have a finite energy in or- 
der to overcome the on-dot Coulomb repulsion and tunnel 
into the dot. Therefore, the conductance of the device is 
suppressed at low temperatures (Coulomb blockade phe- 
nomenon p]). The exceptions are the points of charge 
degeneracy. At these points, two charge states of the 
dot have the same energy, and an electron can hop on 
and off the dot without paying an energy penalty. This 
results in a periodic peak structure in the dependence 
of the conductance G on V g . Away from the peaks, in 
the Coulomb blockade valleys, the charge fluctuations are 
negligible, and the number of electrons N is integer. 

Every time N is tuned to an odd integer, the dot must 
carry a half-integer spin. In the simplest case, the spin is 
S = 1/2, and is due to a single electron residing on the 
last occupied discrete level of the dot. Thus, the quan- 
tum dot behaves as S — 1/2 magnetic impurity imbedded 
into a tunneling barrier between two massive conductors. 
It is known Q since mid-60's that the presence of such 
impurities leads to zero-bias anomalies in tunneling con- 
ductance ||, which are adequately explained M in the 
context of the Kondo effect || . The advantage of the new 
experiments |6| is in full control over the "magnetic impu- 



rity" responsible for the effect. For example, by varying 
the gate voltage, A" can be changed. Kondo effect re- 
sults in the increased low-temperature conductance only 
in the odd- A" valleys. The even- A' valleys nominally 
correspond to the S = spin state (non-magnetic impu- 
rity), and the conductance decreases with lowering the 
temperature. 

Unlike the real atoms, the energy separation between 
the discrete states in a quantum dot is fairly small. 
Therefore, the S — state of a dot with even number 
of electrons is much less robust than the corresponding 
ground state of a real atom. Application of a magnetic 
field in a few-Tesla range may result in a transition to a 
higher-spin state. In such a transition, one of the elec- 
trons residing on the last doubly-occupied level is pro- 
moted to the next (empty) orbital state. The increase in 
the orbital energy accompanying the transition is com- 
pensated by the decrease of Zeeman and exchange ener- 
gies. At the transition point, the ground state of the dot 
is degenerate. Electron tunneling between the dot and 
leads results in mixing of the components of the ground 
state. Remarkably, the mixing involves spin as well as 
orbital degrees of freedom. In this paper we demon- 
strate that the mixing yields an enhancement of the 
low-temperature conductance through the dot. This en- 
hancement can be viewed as the magnetic-field-induccd 
Kondo effect. 

We present the model and theory of electron trans- 
port in the conditions of the field-induced Kondo effect 
in Section [n|. The experimental manifestations of the 
transition observed on GaAs vertical quantum dots and 



carbon nanotubes are described in Section III 
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II. THE MODEL 



E\s,s*) - E\o,o) = KqS - E Z S Z , 



(3) 



We will be considering a confined electron system 
which does not have special symmetries, and therefore 
the single-particle levels in it are non-degenerate. In ad- 
dition, we assume the electron-electron interaction to be 
relatively weak (the gas parameter r s < 1). Therefore, 
discussing the ground state, we concentrate on the tran- 
sitions which involve only the lowest-spin states. In the 
case of even number of electrons, these are states with 

5 = or S = 1. At a sufficiently large level spacing 

6 = e + i — e_i between the last occupied (—1) and the 
first empty orbital level (+1), the ground state is a singlet 
at B = 0. Finite magnetic field affects the orbital ener- 
gies; if it reduces the difference between the energies of 
the said orbital levels, a transition to a state with S = 1 
may occur, see Fig. |l|. Such a transition involves rear- 
rangement of two electrons between the levels n = ±1. 
Out of the six states involved, three belong to a triplet 
5 = 1, and three others are singlets (S = 0). The de- 
generacy of the triplet states is removed only by Zeeman 
energy. The singlet states, in general, are not degenerate 
with each other. To describe the transition between a 
singlet and the triplet in the ground state, it is sufficient 
to consider the following Hamiltonian: 

H dot = tnJ ns d ns - E S S 2 - E Z S Z +E C (N- Nf . 

ns 

(i) 

Here, N = J2 S n 4As i s tne total number of elec- 
trons occupying the levels n = ±1, operator S = 
Unss' d-ns (css'/2) d ns > is the corresponding total spin 
(<x are the Pauli matrices), and the parameters Eg, 
Ez = gfisB, and Ec are the exchange, Zeeman, and 
charging energies respectively [Q. We restrict our at- 
tention to the very middle of a Coulomb blockade val- 
ley with an even number of electrons in the dot (that is 
modelled by setting the dimensionless gate voltage AT to 
Af = 2). We assume that the level spacing 6 is tunable, 
e.g., by means of a magnetic field B: 6 = 5(B), and that 
(5(0) > 2Es (which ensures that the dot is non-magnetic 
for B = 0). 

The lowest-energy singlet state and the three compo- 
nents of the competing triplet state can be labeled as 
\S, S z ) in terms of the total spin S and its z-projection 

S z , 

|i,i) = 4 1T dLi T |o), 

|l,-l) = < u d t _ u |0), (2) 

|i,o) = (< 1T dL u + 4 u dli T ) |o>, 
|o,o) = dl lt dL u |o>, 

where |0) is the state with the two levels empty. Accord- 
ing to ([!]) , the energies of these states satisfy 



where Kq = S — 2E$. Since S > 2E$, the ground state 
of the dot at B — is a singlet |0, 0). Finite field shifts 
the singlet and triplet states due to the orbital effect, 
and also leads to Zeeman splitting of the components of 
the triplet. As B is varied, the level crossings occur (see 
Fig. [I]). The first such crossing takes place at B = B* , 
satisfying the equation 



5{B*)-E Z {B*) = 2E S . 



(4) 



At this point, the two states, |0,0) and |1,1), form a 
doubly degenerate ground state, see Fig. ||. 
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FIG. 1. Typical picture of the singlet-triplet transition in 
the ground state of a quantum dot. 

If leads are attached to the dot, the dot-lead tunneling 
results in the hybridization of the degenerate (singlet and 
triplet) states. The characteristic energy scale To associ- 
ated with the hybridization can be in different relations 
with the Zeeman splitting at field B = B*. 

If Ez(B*) <C To, then the Zeeman splitting between 
the triplet states can be neglected, and at the B = B* 
point all four states @ can be considered as degenerate. 
Theory for this case is presented below in Section [I C . 



This limit adequately describes a quantum dot formed 
in a two-dimensional electron gas (2DEG) at the GaAs- 
AlG aAs i nterface, subject to a magnetic field, see Sec- 
tion III A . Energy Ez can be neglected due to the small- 
ness of the electron g- factor in GaAs. 

Alternatively, the orbital effect of the magnetic field 
(B-dependence of 5) may be very weak due to the re- 
duced dimensionality of the system, while the g-factor 
is not suppressed, yielding an appreciable Zeeman effect 
even in a magnetic field of a mo derate strength. This 
limit of the theory, see Section [ID , corresponds to single- 
wall carbon nanotubes, which have very small widths of 
about 1.4 nm and g — 2.0. Measurements with carbon 



nanotubes are presented in Section [II B . 

In order to study the transport problem, we need to in- 
troduce into the model the Hamiltonian of the leads and 
a term that describes the tunneling. We choose them in 
the following form: 
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Hi — 2^ &C^nks C anksi (5) 
anks 

H T = Yl t <xrm>cl nka d n >, + 'H..C. (6) 
ann' ks 

Here a = R, L for the right/left lead, and n = ±1 for the 
two orbitals participating in the singlet-triplet transition; 
k labels states of the continuum spectrum in the leads, 
and s is the spin index. In writing (@)-@, we had in 
mind the vertical dot device, where the potential creating 
lateral confinement of electrons most probably does not 
vary much over the thickness of the dot || . Therefore we 
have assumed that the electron orbital motion perpendic- 
ular to the axis of the device can be characterised by the 
same quantum number n inside the dot and in the leads. 
Presence of two orbital channels n — ±1 is important for 
the description of the Kondo effect at the singlet-triplet 
transition, that is, when the orbital effect of the magnetic 
field dominates. In the opposite case of large Zeeman 
splitting, the problem is reduced straightforwadly to the 
single-channel one, as we will see in Section II D below. 



A. Effective Hamiltonian 

We will demonstrate the derivation of the effective low- 
energy Hamiltonian under the simplifying assumption M , 



tann' — t a S n 



(7) 



This assumption, on one hand, greatly simplifies the cal- 
culations, and, on the other hand, is still general enough 
to capture the most important physical properties p| . 

It is convenient to begin the derivation by performing 
a rotation II4] in the R-L space 



l^nks 
4*nks 



1 



f tR 




( CRnks \ 




tR J 


V CLnks J 



after which the cj) field decouples: 



Ht = V*L+*?iX>Ld™ + H.c. 



(8) 



(9) 



n ks 



The differential conductance at zero bias G can be re- 
lated, using Eq. (@), to the amplitudes of scattering 
A ns ^ n i s i of the ^-particles 



G = lim dl/dV = - [ -|^% 
v^o ' h \t 2 r + tl 



^ \^-ns—*n' 



(10) 



The next step is to integrate out the virtual transi- 
tions to the states with Af ± 1 electrons by means of the 
Schrieffcr- Wolff transformation or, equivalently, by the 



Brillouin-Wigner perturbation theory. This procedure 
results in the effective low-energy Hamiltonian in which 
the transitions between the states (||) are described by of 
the operators 

ss' 

where V = s z S Z )(S> $ z \ * s the projection oper- 
ator onto the system of states (||). The operators S nn i 
may be conveniently written in terms of two fictitious 
1/2-spins Si : 2- The idea of mapping comes from the one- 
to-one correspondence between the set of states (|j) and 
the states of a two-spin system: 



11,1) 



T1T2 



H,-l> 



I1I2 



|i,o)^=>-^(|Tii a > + Uit2», 
|o,o>^-^(| Ti| 2 )-UiT2». 

We found the following relations: 



S rm — ^ (^1 + — 



1 



(Si 



1 



^mS_ ni „ = V2[Si X S 2 ] = V2T 



(11) 



In terms of Si^, the effective Hamiltonian takes the form: 
H = Yj tA^nks + K (Si ■ S 2 ) - E Z S Z + + H n, (12) 

nks n 

Vn Pnn (Si • S 2 ) (13) 

S_) + 2m (s_„ >n -T)]. 



H n — J (s nn ■ S+) 



Here we introduced the particle and spin densities in the 
continuum: 



Pnn — ^nks^nk'sj s rm' — "^nks ' 



The bare values of the coupling constants are 

J = I = 2V = 2(t 2 L +t R )/E c . (14) 

Note that the Schrieffer- Wolff transformation also pro- 
duces a small correction to the energy gap A between 
the states |1,1) and |0, 0), 



(15) 



so that K differs from its bare value Kq, see ([|). However, 
this difference is not important, since it only affects the 
value of the control parameter at which the singlet-triplet 
transition occurs, but not the nature of the transition. 
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We did not include into (pL 2|) - ([1 3|) the free-electron 
Hamiltonian of the ^-particles [see Eq. (13)], as well as 
some other terms, that are irrelevant for the low energy 
renormalization. The contribution of these terms to the 
conductance is featureless at the energy scale of the order 
of To (see the next section), where the Kondo resonance 
develops. 

At this point, it is necessary to discuss some approxi- 
mations tacitly made in the derivation of (|l2])-(|l5|). First 
of all, we entirely ignored the presence of many energy 
levels in the dot, and took into account the low-energy 
multiplet (]|) only. The multi-level structure of the dot is 
important at the energies above S, while the Kondo effect 
physics emerges at the energy scale well below the single- 
particle level spacing |l3|]. The high-energy states result 
merely in a renormalization of the parameters of the effec- 
tive low-energy Hamiltonian. One only needs to consider 
this renormalization for deriving the relation between the 
parameters ti and of the low-energy Hamiltonian (Q) , 
(||) and (||) and the "bare" constants of the model defined 
in a wide bandwidth ep. On the other hand, using the 
effective low-energy Hamiltonian, one can calculate, in 
principle, the observable quantities such as conductance 
G(T) and other susceptibilities of the system at low tem- 
peratures (T <C 6), and establish the relations between 
them, which is our main goal. 

Note that the Hamiltonian (p"2|)-(p"3|) resembles that 
of the two-impurity Kondo model, for which H n — 
J n (s nn ■ S+) +/ (s_ ni „ • S_) and the parameter K char- 
acterizes the strength of the RKKY interaction 0|. It 
is known that the two-impurity Kondo model may un- 
dergo a phase transition at some special value of K [0 . 
At this point, the system may exhibit non- Fermi liquid 
properties. However, one can show [Io| , using general 
arguments put forward in (l4|, that the model (|l^)-(|l3|) 
does not have the symmetry that warrants the existence 
of the non Fermi liquid state. This allows one to apply 
the local Fermi liquid description [jl5| to study the prop- 
erties of the system at T = 0. In the next section, we will 
concentrate on the experimentally relevant perturbative 
regime. 



B. Scaling Analysis 



To calculate the differential conductance in the leading 
logarithmic approximation, we apply the "poor man's" 
scaling technique The procedure consists of a per- 
turbative elimination of the high-energy degrees of free- 
dom and yields the and yields the set of scaling equations 



dJ/dC = v ( J 2 + I 2 ) , 
dl/dC = 2vI(J + V), 
dV/dC = 2vl 2 



(16) 



for the renormalization of the coupling constants with the 
decrease of the high energy cutoff D. Here C = \n(S/D), 
and v is the density of states in the leads; the initial value 
of D is D — S, see the discussion after Eq. (|l5|). The ini- 
tial conditions for (|l6|), J(0), 1(0), and V(0) are given 
by Eq. (|l~4|). The scaling procedure also generates non- 
logarithmic corrections to K . In the following we absorb 
these corrections in the re-defined value of K . Equations 
( pf ) are valid in the perturbative regime and as long as 

£>» \K\,E Z ,T. 

At certain value of C, C = Co = ln(<5/T ), the inverse 
coupling constants simultaneously reach zero: 

l/J(£o) = l//(£o) = l/7(£o) = 0. 

This defines the characteristic energy scale of the prob- 
lem: 



T = S exp [-Tq/vJ] 



(17) 



Here To is a parameter that depends on the initial con- 
ditions and should be found numerically. We obtained 
t = 0.36 (see Fig. |). 



x(x) y(i) 




FIG. 2. Numerical solution of the scaling equations. 
The RG equations @ are rewritten in terms of the 
new variable t = v J(0) \n(S/D) and the new functions 
x(r) = J(0)/J(t), y(r) = /(0)/J(r), z{r) = V(0)/V(r) 
as dx/dr — —(1 + x 2 /y 2 ), dy/dr = —(2y/x + y/z), 
dz/dr = —Az 2 /y 2 . The three functions reach zero simulta- 
neously at t — to = 0.36. 

It is instructive to compare To with the Kondo tem- 
perature T£- dd in the adjacent Coulomb blockade valleys 
with N = odd. In this case, only electrons from one 
of the two orbitals n = ±1 are involved in the effec- 
tive Hamiltonian, which takes the form of the 1-channel 
S = 1/2 Kondo model with the exchange amplitude 
Jodd = 4(i| + t%)/E c = 2J, see Eq. (fil). Therefore, 



Ddd 
T^odd 
1 K 



is given by the same expression (|17|) as To, but 
with To = 1/2. For realistic values of the parameters 
T = 300 mK, S = 3 meV we obtain T° dd rs 120 mK. 
This estimate is in a reasonable agreement with the ex- 
perimental data, see Section III A below. 

The solution of the RG equations ( |l6[ ) can now be ex- 
panded near C = Co- To the first order in Co — C = 
In D/ To, we obtain 



1 



Va a- i 



vJ(C) vl(£) 2vV{C) 



= (A + l)ln(D/T ), (18) 
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where 



A = 2 + \/5«4.2. 



It should be emphasized that, unlike To, the constant A 
is universal in the sense that its value is not affected if 
the restriction (0) is lifted Q . 

Eq. ( |l8| ) can be used to calculate the differential con- 
ductance at high temperature T 3> \K\ ,Ez,T n . In this 
regime, the coupling constants are still small, and the 
conductance is obtained by applying a perturbation the- 
ory to the Hamiltonian (|l2|)-(|l3|) with renormalized pa- 
rameters (|i"8|), taken at D = T, and using (|l0|). This 
yields 



G/Gq — 



A 



MT/T )] 



2 ' 



(19) 



where 



A= (3^ 2 /8) (A + ir 2 1 + A+(A-1) 



0.9 



is a numerical constant, and 



4e 2 



2t T .t 



L L R 



j.2 I j.2 
L L ^ L R 



(20) 



As temperature is lowered, the scaling trajectory Jl6| ) 
terminates either at D ~ max{|A|,Fz} 3> To, or when 
the system approaches the strong coupling regime D ~ 
To ^> \K\,Ez- It turns out that the two limits of the 
theory, Ez <C T and Ez ^S> To, describe two distinct 
physical situations, which we will discuss separately. 



C. Singlet-triplet transition 

In this section, we assume that the Zeeman energy is 
negligibly small compared to all other energy scales. At 
high temperature T 3> |if|,To, the conductance is given 
by Eq. (|l9|). At low temperature T <\K\ and away from 
the singlet-triplet degeneracy point, \K\ ^> To, the RG 
flow yielding Eq. (|l9| ) terminates at energy D ~ \K\. On 
the triplet side of the transition (K <C —To), the two 
spins Si,2 are locked into a triplet state. The system is 
described by the effective 2-channel Kondo model with 
5 = 1 impurity, obtained from Eqs. (p^|)-(|l3|) by pro- 
jecting out the singlet state and dropping the no longer 
relevant potential scattering term: 



^triplet - ^2£ki> nks 



4>n 



(21) 



here J is given by the solution J(C) of Eq. (|l6|), taken at 
£ = £* = ln{S/\K\), which corresponds to D = \K\. 

As D is lowered below \K\, the renormalization of the 
exchange amplitude J is governed by the standard RG 
equation |16 



dJ/dC = vJ 1 

where C = \jx{5/D) > C*. Eq. 
with the result 



(22) 

2|) is easily integrated 



l/vJ{C) - l/v,]{C*) =£-£*. 

This can be also expressed in terms of the running band- 
width D and the Kondo temperature 

T k = |K|exp[-EW(£*)] 

as l/vJ{£) = \n{D/T k ). 

Obviously, T k depends on \K\. Using asymptotes of 
J{C), see Eq. @, we obtain the scaling relation 



T k /T = (T /\K\Y 



(23) 



Eq. ( p3| ) is valid not too far from the transition point, 
where the inequality 



1 < |A|/T « (S/ToT, Ai«0.24 



(24) 



is satisfied. Here, \x is a numerical constant, which de- 
pends on To, and therefore is not universal [see the remark 
after Eq. (|l8|)]. For larger values of \K\ (but still smaller 
than (5), T fe oc 1/\K\ Q. Finally, for \K\ = S, T k is given 
by Eq. (|l7|) with t = 1 . According to (^3|) , T k decreases 
very rapidly with \K\. For example, for To = 300 mK 
and 6 = 3 meV Eq. (g3|) describes fall of T k by an order 
of magnitude within the limits of its validity (|24|). For 
| if | = 6 one obtains T k as 5 mK, which is well beyond 
the reach of the present day experiments. 

For a given \K\, T < \K\ < 8, the differential conduc- 
tance can be cast into the scaling form, 



G/G =F(T/T k ) 



(25) 



where F (x) is a smooth function that interpolates be- 
tween F(0) = 1 and F (x > 1) = (tt 2 /2) (lna;)" 2 . It 
coincides with the scaled resistivity 

F(T/T K ) = p(T/T K )/p(0) 

for the symmetric two-channel S = 1 Kondo model. The 
conductance at T = (the unitary limit value), Go, is 
given above in Eq. (po|) . 



FIG 



-r ° r ^ 

3. Linear conductance near a singlet-triplet transi- 
tion. At high temperature G exhibits a peak near the transi- 
tion point. At low temperature G reaches the unitary limit at 
the triplet side of the transition, and decreases monotonously 
at the singlet side. The two asymptotes merge at K ^> T, T . 
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On the singlet side of the transition, K Tq, the scal- 
ing terminates at D ~ K, and the low-energy effective 
Hamiltonian is 



^singlet = Yl/^nks^nks ~ J V Y1 



nks 



where V is V(£) [see Eq. (|l|)] taken at C = C* . The 
temperature dependence of the conductance saturates at 
T <C K, reaching the value 



G/Gq 



B 



[HK/T )Y 



B 



3tt A - 1 



0.5. (26) 



Note that at T = Eqs. fl25J) and predict different 
dependence on the parameter K which is used for tuning 
thorugh the transition. At positive K , conductance de- 
creases with the increase of K; at K -C — Tq conductance 
G = Go and does not depend of K. Although there is no 
reason for the function G(K) to be discontinious |l0[, it 
is obviously a non- analytical function of K, see Fig. |3|. 

The above results are for the linear conductance G. At 
T = 0, G is a monotonous function of K, at high tem- 
perature T 3> To the conductance develops a peak at the 
singlet-triplet transition point K = 0. We now discuss 
shortly out-of-equilibrium properties. When the system 
is tuned to the transition point K — 0, the differential 
conductance dl/dV exhibits a peak at zero bias, whose 
width is of the order of To. For finite K the peak splits 
in two, located at finite bias eV = ±K. The mecha- 
nism of this effect is completely analogous to the Zee- 
man splitting of the usual Kondo resonance Q, \K\ 
is the energy cost of the processes involving a singlet- 
triplet transition. This cost can be covered by applying 
a finite voltage eV = ±K, so that the tunneling electron 
has just the right amount of extra energy to activate the 
singlet-triplet transition prosesses described by the last 
two terms in ([l3]) . The split peaks gradually disappear at 
large \K\ due to the nonequilibrium- induced decoherence 

in, H. 



D. Transition Driven by Zeeman Splitting 



If the Zeeman energy is large, the RG flow (16) ter- 
minates at D ~ Ez- The effective Hamiltonian, valid at 
the energies D < Ez is obtained by projecting (p^)-(p^) 
onto the states |1,1) and |0, 0). These states differ by 
a flip of a spin of a single electron (see Fig. ||), and are 
the counterparts of the spin-up and spin-down states of 
5 = 1/2 impurity in the conventional Kondo problem. It 
is therefore convenient to switch to the notations 



|i,i) = |T>, |o,o) = ||), 



(27) 



and to describe the transitions between the two states in 
terms of the spin-like operator 



built from the states (E 



4- 



[l,l>=it) [o,o> = |*> 

FIG. 4. The ground state doublet in case of a large Zee- 
man splitting. The states |1, 1} and |0, 0) differ by flipping a 
spin of a single electron (marked by circles). 

Projecting onto the sates (p7|), we obtain from (|l2|)- 

nks 

+ E [ J <n (S z + 1/2) + Vnp nn (> - 1/4)] (28) 

n 

—I (sf_ 1 S~ + H.c. 



where A was introduced above in Eq. (|15|). It is now 
convenient to transform (Pq) to a form which is diago- 
nal in the orbital indexes n. This is achieved simply by 
relabeling the fields according to 



1P-l,k,t = b k;h V'+l.fc.l = 



(29) 



which yields 
H = H + AS Z 



+v a s z a + J z s z a s z + -J A (,,; s- 

+V b s z b + J' 2 s z b S z , 



"5+ 



(30) 



where Hq is a free-particle Hamiltonian for a, b elec- 
trons, and s a is the spin density for a electrons, s a — 
Sfcfc'ss' a L (°w/2) a k's' (with a similar definition for 
s b ). The coupling constants in (30), 



V a = (J - V)/2, J Z = J 
V b = (J+V)/2, J' Z = J- 



- 2V, J± 
2V 



2L 



(31) 



are expressed through the solutions of the RG equations 
@ taken at C = £** = \n{5/E z ). 

The operators in 6-dependent part of (^) are not rel- 
evant for the low energy renormalization. At low enough 
temperature (satisfying the condition \n{T/Tz) <C 
(^J^) _1 , (^Vfc) -1 , where Tz is the Kondo temperature), 
their contribution to the conductance becomes negligi- 
ble compared to the contribution from the a-dependent 
terms. This allows us to drop the 6-dependent part of 
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(pO|). Suppressing the (now redundant) subscript of the 
operators s l a , we are left with the Hamiltonian of a one- 
channel S = 1/2 anisotropic Kondo model, 

Jj 



H = H + AS Z + V a s z + J z s z S z 



-(s+S- +H.c). 



(32) 



Eq. ( p2] ) emerged as a limiting case of a more general 
two-channel model (^2|)-(^3[). It should be noticed, how- 
ever, that the same effective Hamiltonian (|3^) appears 
when one starts with the single-channel model from the 
very beginning Jn|. 

A finite magnetic field singles out the z-direction, so 



that the spin-rotational symmetry is absent in (|12|)-(13). 
This property is preserved in (p2[). Indeed, even for 



J z = J±, Eq. ( p2|) contains term V,ps z which has the 
meaning of a magnetic field acting locally on the conduc- 
tion electrons at the impuirity site. The main effect of 
this term is to produce a correction to A, through creat- 
ing a non-zero expectation value (s z ) |po[ . This results 
in a correction to A. Fortunately, this correction is not 
important, since it merely shifts the degeneracy point. 
In addition, this term leads to insignificant corrections 
to the density of states . 

Let us now examine the relation between J z and J±. 
It follows from Eqs. (|l4|) and @, that at E z = 6 the ex- 
change is isotropic: J z = J±. Moreover, it turns out that 
if Ez is so close to 5, that Eqs. ( |l6| ) can be linearized 
near the weak coupling fixed point C = 0, the correc- 
tions to J z , Jj_ are such that the isotropy of exchange is 
preserved: 



Jz = J± 



2J(0) 



1 



-vJ(0)hx(S/E z ) 



(33) 



where J(0) is given by (|lj). This expression is valid 
as long as the logarithmic term in the r.h.s. is small: 
vJ(0)ln(6/E z ) < 1. Using (§|) and (0), one obtains 
the Kondo temperature Tz, which for J z = J± is given 
by 

T z = Ezexp[-l/vJ z ] = E z (5/Ez) 3/s (To/S) 1/2t ° 

Note that for E z = 6, Tz coincides with the Kondo tem- 
perature Tg- dd in the adjacent Coulomb blockade valleys 
with odd number of electrons [|l9| , see the discussion after 
Eq. @ above. 

Note that the anisotropy of the exchange merely af- 
fects the value of Tz (which can be written explicitely 
for arbitrary J z and J± |2l[| ). In the universal regime 
(when T approaches Tz), the exchange can be consid- 
ered isotropic. This is evident from the scaling equations 

fi 

dJ z /dC = vj\, dJ ± /dC = vJ z J±, £> £** (34) 

where C > £** = \n(S/Ez), whose solution approaches 
the line J z — J \ at large C 



According to the discussion above, the term V^s z in 
Eq. ( p2| ) can be neglected. As a results, (32) acquires the 
form of the anisotropic Kondo model, with A playing 
the part of the Zeeman splitting of the impurity levels. 
This allows us to write down the expression for the linear 
conductance at once. Regardless the initial anisotropy of 
the exchange constants in Eq. |3^), the conductance for 
A = in the universal regime (when T approaches Tz or 
lower) is given by 



G = G 0Z f(T/T z ) 



(35) 



where f(x) is a smooth function interpolating between 
/(0) = 1 and f(x > 1) = (37r 2 /16)(lnx)- 2 . Function 
/ (T/Tz) coincides with the scaled resistivity for the one- 
channel S — 1/2 Kondo model and its detailed shape is 
known from the numerical RG calculations p^]. The 
conductance at T = 0, 



Gqz — —r~ 



2e_ 2 
h 



(36) 



is by a factor of 2 smaller than G [see Eq. (poj)]; Go 
includes contributions from two channels and therefore 
is twice as large as the single-channel result ([56]). At 
finite A ^> Tz, the scaling trajectory ( |34| ) terminates 
at D ~ A. As a result, at T < A the conductance is 
temperature-independent, and for J z = Ji 



G = G 0Z f(A/T z ) = G, 



az 



3tt 2 /16 

MNTz)f 



The effect of the de-tuning of the magnetic field from the 
degeneracy point A / on the differential conductance 
away from equilibrium is similar to the effect the mag- 
netic field has on the usual Kondo resonance Q , [0 . For 
example, consider the case, relevant for the experiments 
on the carbon nanotubes, see section III B , when the ex- 
change energy Es [see Eq. d3)] is negligibly small. When 
sweeping magnetic field from B = — oo to B = +oo, the 
degeneracy between the singlet state of the dot and a 
component of the triplet is reached twice, at B = B* 
and B = —B*, when \E Z \ ~ S. If the field is tuned to 
B = ±B* , then the differential conductance dl/dV has 
a peak at zero bias. At a finite difference \B\ — \B*\ this 
peak splits in two located at eV — ±g/i£(|£>| — |£?*|). 



III. EXPERIMENTS 

A. GaAs quantum dots 

Here, we discuss the case of a quantum dot with 
N = even in a situation where the last two electrons 
occupy either a spin singlet or a spin triplet state. The 
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transition between singlet and triplet state is controlled 
with an external magnetic field. The range of the mag- 
netic field is small (B ~ 0.2 T, g^sB ~ 5 fiV) such 
that the Zeeman energy can be neglected and that the 
triplet state is fully degenerate p3||. The theory for this 



situation was described in section [I C 




FIG. 5. (a) Cross-section of rectangular quantum 

dot. The semiconductor material consists of an undoped 
AlGaAs(7nm)/InGaAs(12nm)/AlGaAs(7nm) double barrier 
structure sandwiched between n-doped GaAs source and 
drain electrodes. A gate electrode surrounds the pillar and 
is used to control the electrostatic confinement in the quan- 
tum dot. A dc bias voltage, V, is applied between source and 
drain and current, I, flows vertically through the pillar. The 
gate voltage, V g , can change the number of confined electrons, 
N, one-by-one. A magnetic field, B, is applied along the ver- 
tical axis. 

(b) Scanning electron micrograph of a quantum dot with di- 
mensions 0.45 x 0.6 /jm 2 and height of ~ 0.5 \im. 

The quantum dot has the external shape of a rect- 
angular pillar (see Fig. 0) and an internal confinement 
potential close to a two-dimensional ellipse || . The tun- 
nel barriers between the quantum dot and the source 
and drain electrodes are thinner than in other devices 
such that higher-order tunneling processes are enhanced. 
Fig. |^ shows the linear response conductance G versus 
gate voltage V g , and magnetic field B. Dark regions 
have low conductance and correspond to the regimes 
of Coulomb blockade for N = 3 to 10. Light stripes 
represent Coulomb peaks as high as ~ e 2 /h. The in- 
dependence of the first two lower stripes reflects the 
ground-state evolution for N = 3 and 4. Their simi- 
lar B-evolution indicates that the 3rd and 4th electron 
occupy the same orbital state with opposite spin, which 
is observed also for N — 1 and 2 (not shown) . This is not 
the case for N = 5 and 6. The N = 5 state has S = 1/2, 
and the corresponding stripe shows a smooth evolution 
with B. Instead, the stripe for N — 6 has a kink at 
B = B* ps 0.22 T. From earlier analyses || and from 
measurements of the excitation spectrum at finite bias 
V this kink is identified with a transition in the ground 
state from a spin-triplet to a spin-singlet. 

Strikingly, at the triplet-singlet transition (see Fig. |^) 
we observe a strong enhancement of the conductance. In 
fact, over a narrow range around 0.22 T, the Coulomb 
gap for TV = 6 has disappeared completely. Note that 
the change in greyscale along the dashed line in Fig. ^ 
represents the variation of the conductance with the tun- 
ing parameter K, see Fig. ||. 



To explore this conductance anomaly, Fig. ^(a) shows 
the differential conductance, dl/dV versus V, taken at 
B and V g corresponding to the intersection of the dotted 
line and the bright stripe (B = B*) in Fig. ^. The height 
of the zero-bias resonance decreases logarithmically with 
T [see Fig. 0(b)]. These are typical fingerprints of the 
Kondo effect. From FWHMw 30 fiV « k B T , we esti- 
mate T m 350 mK. Note that k B T /g^ B B* m 6 so 
that the triplet state is indeed three-fold degenerate on 
the energy scale of To; this justifies an assumption made 
in Section fl C| above. Also note that some of the traces 
in Fig. ^(a) show small short-period modulations which 
disappear above ~ 200 mK. These are due to a weak 
charging effect in the GaAs pillar above the dot [E4j . 




Magnetic field (T) 



FIG. 6. Gray-scale representation of the linear conduc- 
tance G versus the gate voltage V g and the magnetic field B. 



White stripes denote conductance peaks of height 



Dark regions of low conductance indicate Coulomb blockade. 
The N = 6 ground state undergoes a triplet-to-singlet tran- 
sition at B — B* ~ 0.22 T, which results in a conductance 
anomaly inside the corresponding Coulomb gap. 

For N = 6 the anomalous T-dependence is found only 
when the singlet and triplet states are degenerate. Away 
from the degeneracy, the valley conductance increases 
with T due to thermally activated transport. For N = 5 
and 7, zero-bias resonances are clearly observed [see in- 
sets to Fig. 0(a)] which are related to the ordinary spin- 
1/2 Kondo effect. Their height, however, is much smaller 
than for the singlet-triplet Kondo effect. 

We now investigate the effect of lifting the singlet- 
triplet degeneracy by changing B at a fixed V g corre- 
sponding to the dotted line in Fig. ^|. Near the edges of 
this line, i.e. away from B* , the Coulomb gap is well 
developed as denoted by the dark colours. The dl/dV vs 

V traces still exhibit anomalies, however, now at finite 

V [see Fig. |). For B = 0.21 T we observe the singlet- 
triplet Kondo resonance at V = 0. At higher B this 
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resonance splits apart showing two peaks at finite V, in 
agreement with the discussion above (see Section II C ). 
For B w 0.39 T the peaks have evolved into steps which 
may indicate that the spin-coherence associated with the 
Kondo effect has completely vanished. The upper traces 
in Fig. H for B < 0.21 T, also show peak structures, 
although less pronounced. 



0.75 




0.75 



-0.5 0.5 

Bias voltage (mV) 




10 1 10 2 10 J 

Temperature (mK) 



FIG. 7. (a) Kondo resonance at the singlet-triplet transi- 
tion. The dl/dV vs V curves are taken at V g = —0.72 V, 
B = 0.21 T and for T = 14,65,100,200,350,520, and 
810 mK. Kondo resonances for N = 5 (left inset) and N = 7 
(right inset) are much weaker than for N = 6. 
(b) Peak height of zero-bias Kondo resonance vs T as obtained 
from (a). The line demonstrates a logarithmic T-dependence, 
which is characteristic for the Kondo effect. The saturation 
at low T is likely due to electronic noise. 



very weak, because the tube diameter (~ 1.4 nm) is 
an order of magnitude smaller than the magnetic length 
lg = (h/eB) 1 / 2 ~ 10 nm at a typical maximum labo- 
ratory field of 10 T. On the other hand, the ^-factor is 
close to its bare value of g = 2, compared with g = 0.44 
in GaAs. Hence the magnetic response of a nanotube 
dot is determined mainly by Zeeman shifts. As a result, 
the spins of levels in nanotube dots are easily measured 
p6[ , pS|, and the ground state is usually (though not al- 
ways [|26[) found to alternate regularly between an S = 
singlet for even electron number N and an S = 1/2 dou- 
blet for odd N j28|, p9[ |. Moreover, singlet-triplet transi- 
tions in nanotubes are likely to be driven by the Zeeman 
splitting rather than o rbita l shifts, corresponding to the 



theory given in section [II D 




FIG. 9. Schematic of a nanotube quantum dot, incorpo- 
rating an atomic force microscope image of a typical device 
(not the same one measured here.) Bridging the contacts, 
whose separation is 200 nm, is a 2 nm thick bundle of sin- 
gle-walled nanotubes 




Bias voltage (mV) 

FIG. 8. dl/dV vs V characteristics taken along the 
dotted line in Fig. ^ at equally spaced magnetic fields 
B = 0.11, 0.13, 0.39 T. Curves are offset by 0.25 e 2 /h. 

B. Carbon nanotubes 

The situation in quantum dots formed in single-wall 
carbon nanotubes |||, @, Q, Q is rather differ- 
ent from that in semiconductor quantum dots. In nan- 
otubes the effect of magnetic field on orbital motion is 



We discuss here the characteristics of a single-walled 
nanotube device with high contact transparencies, which 
were presented in more details in |29[ . The source and 
drain contacts are gold, evaporated on top of laser 
ablation-grown nanotubes |3(| deposited on silicon diox- 
ide. The conducting silicon substrate acts as the gate, 
as illustrated in Fig. ^|. At room temperature the linear 
conductance G is 1.6 e 2 /h, almost independent of gate 
voltage V g , implying the conductance-dominating nan- 
otube is metallic and defect free, and that the contact 
transmission coefficients are not much less than unity. 
At liquid helium temperatures regular Coulomb block- 
ade oscillations develop, implying the formation of a sin- 
gle quantum dot limiting the conductance. However, the 
conductance in the Coulomb blockade valleys does not 
go to zero, consistent with high transmission coefficients 
and a strong coupling of electron states in the tube with 
the contacts. 

Fig. [l(] shows a grayscale plot of dl/dV versus V and 
V g over a small part of the full V g range at B = 0. A regu- 
lar series of faint "Coulomb diamonds" can be discerned, 
one of which is outlined by white dotted lines. Each di- 
amond is labeled either E or O according to whether N 
is even or odd respectively, as determined from the ef- 
fects of magnetic field. Superimposed on the diamonds 
are horizontal features which can be attributed to higher- 
order tunnelling processes that do not change the charge 
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on the dot and therefore are not sensitive to V g . 




2.35 2.40 2.45 2.50 
KM 



FIG. 10. Grayscale plot of differential conductance dl jdV 
(darker = more positive) against bias V and gate voltage V g at 
a series of magnetic fields and base temperature (~ 75 mK). 
Labels 'E' and 'O' indicate whether the number of electrons 
TV in the dot is even or odd (see text). 



In Fig. pj](a) we concentrate on an adjacent pair of E 
and O diamonds in a magnetic field applied perpendic- 
ular to the tube. At 73 = the diamond marked with 
an 'O' has narrow ridge of enhanced dl/dV spanning it 
at V = 0, while that marked with an 'E' does not. An 
appearence of a ridge at zero bias is consistent with for- 
mation of a Kondo resonance which occurs when N is odd 
(O) but not when it is even (E). This explanation is sup- 
ported by the logarithmic temperature dependence of the 
linear conductance in the center of the ridge, as indicated 
in the inset to Fig. |l^(b). At finite 73, each zero- bias 
ridge splits into features at approximately V = ±Ez/e 
as expected for Kondo resonances B . 




(b) 

















2.35 T 


V = -0.322 V 
i . i 





-0.2 0.0 
V (mV) 



0.2 



FIG. 11. (a) Evolution with magnetic field of adjacent 
even (E) and odd (O) features of the type seen in Fig. Fuj. 
(b) dl/dV vs V traces at the center of the E region, at 
V a = -0.322 V. The trace at 73 = 73* = 1.18 T (bold line) 
corresponds to the dotted line in (a). The traces are offset 
from each other by 0.4 e 2 /h for clarity. 



On the other hand, in the E diamond the horizon- 
tal features appear at a finite bias at B = 0. The ori- 
gin of these features can be infered from their evolution 
with a magnetic field: while the ridge in the O region 
splits as 73 increases, the edges of the E 'bubble' move 
towards V = 0, finally merging into a single ridge at 
B = B* = 1.18 T. Fig. [U](b) shows the evolution with 
B of the dl/dV vs V traces from the center of the E re- 
gion, and the appearance of a zero-bias peak at around 
B = 1.18 T (bold trace). This matches what is expected 
for a Zeeman-driven singlet-triplet transition in the N = 
even dots (Section |II Ej ). Further evidence that the peak 
is a Kondo resonance is provided by its temperature de- 
pendence [Fig. |l2|(a)] , which shows an approximately log- 
arithmic decrease of the peak height (the linear conduc- 
tance) with T shown in Fig. |l2|(b). 




-0.2 0.0 
V (mV) 



0.6 



T (K) 



FIG. 12. (a) Temperature dependence at 73 = 73*. Here 
T= base (bold line), 100, 115, 130, 180, 230, and 350 mK. 
(b) Temperature dependence of the linear conductance G 
(dl/dV at V = 0) at 73 = 73* = 1.18 T. For comparison, 
G{T) in the center of one of the O-type ridges at 73 = is 
shown in the inset. 



Based on this interpretation we can deduce that for 
this particular value of N the energy gap separating the 
singlet ground state and the lowest-energy triplet state 
is Ao = g^sB* w 137 fieV. At other even values of 

the lowest visible excitations range in energy up to 
~ 400 [ieV . For this device E c ~ 500 [ieV . The energy 
gaps are therefore comparable with the expected single- 
particle level spacing 5, which is roughly equal to Eq/3 
in a nanotube dot |27l] . 



Note that the ridges at finite bias in Fig [Tl] in E valley 
are more visible at 73 = than at B — 0.59 T, halfway 
towards the degeneracy point. A possible explanation is 
that at 73 = the triplet is not split, and all its com- 
ponents should be taken into account when calculating 
dl/dV at 73 = 0. This results in an enhancement of 
dl/dV at 73 = 0, eV — S, as compared to the value ex- 
pected from the effective model of Section 
valid in the vicinity of 73 = 73*. 



[ID 
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CONCLUSION 



Even a moderate magnetic field applied to a quantum 
dot or a segment of a nanotube can force a transition 
from the zero-spin ground state (S — 0) to a higher-spin 
state (S = 1 in our case). Therefore, the magnetic field 
may induce the Kondo effect in such a system. This is 
in contrast with the intuition developed on the conven- 
tional Kondo effect, which is destroyed by the applied 
magnetic field. In this paper we have reviewed the ex- 
perimental and theoretical aspects of the recently studied 
magnetic-field-induced Kondo effect in quantum dots. 
Clearly there is more territory to be explored in the re- 
markably tuneable systems. 
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